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ABSTRACT 



Context. It is usual in helioseismology to remove unwanted instrumental low-frequency trends by applying high-pass filters to the 
time series. However, the choice of the filter is very important because it can keep the periodic signals throughout the spectrum. At the 
same time, these filters should not introduce any spurious elfects on the remaining signal, which would modify the periodic signatures 
in the Fourier domain. 

Aims. One of the most used filters is the so-called backwards dilference filter that can be applied when the time series are regularly 
sampled. The problem of this filter is that the amplitudes and the phases of the periodic signals in the Fourier domain are modified and, 
therefore, their extracted amplitudes are biased when the frequency of the periodic signals decreases. The objective of this research 
note is to give a correction that could be applied to the power spectrum to over come this problem. 

Methods. We properly derive the transfer function of the backwards dilference filter and we show how this can be applied to the 
power spectrum to correct the amplitudes. For the amplitude spectrum and the time series, we also derive a correction for the phase. 
Results. The amplitudes of the periodic signals in the resultant power spectrum density are corrected from the cut-olf frequency down 
to zero. 

Key words. Methods: data analysis - Methods: filters - Stars: oscillations ~ Sun: helioseismology 



1. Introduction 

One of the most common problems to be faced when analyz- 
ing oscillatory components of an observed signal is the presence 
of low-frequency drifts which can be of instrumental or natural 
origin. These low-frequency trends introduce a background on 
the Fourier domain which reduces the signal-to-noise ratio (S/N) 
of searched periodic signals. To avoid this, high-pass filters are 
commonly used. The ideal filter is one that reduces the effect of 
the trend while preserving the periodic signals even at low fre- 
quencies. Moreover, it is important to keep all the parameters of 
periodic signals untouched, including their amplitudes. 

In this research note we study a very easy-to-use filter: the 
backwards difference filter w hich is commonly used in he lio- 
seismology (see for example IChristensen-DalsgaardI l2002l for 
a complete review on this research topic) to reduce the low- 
frequency unwanted stochastic contribution of regularly sampled 
observations. We also deduce the transfer function of this fil- 
ter. This transfer function allows us to correct the filtered power 
spectrum in order to retrieve the correct amplitudes of the peri- 
odic signals in the spectrum. Finally, we determine a correction 
for the resultant phase for when we are interested in the ampli- 
tude spectrum or in the time series. 



2. Backwards difference filter 

2.1. Definition 

The backwards difference filter (hereafter called BDf) can be ob- 
tained by substituting the value of every data point of a regularly 
sampled time series by the difference of two consecutive points. 
Let {/„) be a set of regularly sampled measurements of a func- 
tion /(f), such that /„ = f{t„) with t„ - fo + nAt, where At is the 
sampling rate. By definition, the filtered data set {6f„] is such 
that 

Sfri - fn+\ - fn (1) 

which can be seen as the discretization of the function 

6f(t)^f(t + At)-f(t) (2) 

Thus, for a given Af, the BDf can be defined as a linear operator, 
6, that can be applied to any function /(f) using Eq.|2] 

2.2. Transfer function of the backwards difference filter 

Denoting /(w) - !F[/(f)] the Fourier transform of the function 
/(f), we can write that: 

6f(cj) = r[/(f + Af)-/(f)] 
- [exp(/dLtAf) - 1] f(cD) 
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/ wAf \ ~ I LoAt + 7T 
2 sin I — I /(w) exp I i ^ 



(3) 
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We deduce that the power spectrum of the original function /(f) 
and the filtered one, 6f(t), are linked by: 



(4) 



where Q((l>) is the transfer function of the BDf. It can be defined 
from the previous expressions by: 
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(5) 



where we have introduced the cut-off frequency ojc - n/At. 
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Fig. 1. Transfer function Q{a>) (defined in Eq. |5]l of the back- 
wards difference filter. The frequency w has been normalized to 
the cut-off frequency ojc - n/At. 

This transfer function shows the effect of the filter on the 
amplitude of the power spectrum (see Fig. [1]). Indeed, at lower 
frequencies, the amplitudes tend towards zero. To recover the 
correct amplitudes in the filtered power spectrum, we need to 
divide it by Qicj): 



\r[f(tw = Qio^r \r[6f(t)]f 



3. Example 



(6) 



To illustrate the behavior of the BDf and to show the advantage 
of such a filter, we have simulated a simple observation contain- 
ing low-frequency drifts and periodic signals (see Fig. |2]i. The 
length of the simulated series is T = 150000 s with a constant 
sampling rate. At = Is. We have simulated three sine waves with 
periods of 1 000, 5 000 and 10000 s, corresponding to frequen- 
cies of 10"^, 2x10""* and 10""* Hz, in a way to cover a large range 
of frequencies. Their amplitudes are all equal to 1 000 arbitrary 
units. A low-frequency trend has been added to mimic the be- 
havior of real observations. Thus, we have used a 7*-order poly- 
nomial function fitted to the light curve of the Procyon data re- 
covered by the MOSlQ satellite. This data set ca n be found in the 
MOST public data archiv^ and it is described in lMatthews et alj 

JIool . 



' A Canadian Space Agency mission, jointly operated by Dynacon 
Inc., the University of Toronto Institute for Aerospace Studies and the 
University of British Columbia, with the assistance of the University of 
Vienna 

^ http://www.astro.ubc.ca/MOST/data/data.html 
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Fig. 2. Simulated time series of the sum of three sine waves and a 
background that mimics a standard asteroseismic observed back- 
ground. 



The power spec tral density (PSD) - calibrated as a single- 
sided spectrum, see iPress et al] (1 1 9921) - of the simulated time 
series is plotted in Fig. [3^. The low-frequency trend introduces 
a background in the power spectrum reducing the S/N of the 
sine waves as the frequency decreases when it is compared to 
the PSD of the three sine waves computed without any trend 
(cf. Fig.[3}3). The application of the BDf reduces the effect of the 
unwanted background (cf. Fig. [3]:). Indeed the S/N ratio increase 
compared to the raw PSD. However, the amplitudes of the sine 
waves are modified by the transfer function of the filter. Thus, 
when we divide this PSD by the Q(cj) factor, we retrieve the 
correct amplitudes (cf. Fig.[3}l)- 

Another way sometimes used to filter slow trends is to re- 
move a smoothing of time series. To make a comparison, we 
have also plotted the PSD of our artificial data after applying 
such a filter In this example, we have smoothed our time series 
by a boxcar function with a width of T/ 10 and subtracted it from 
the data. Figure[3^ clearly shows that, in this last case, the back- 
ground of the PSD is dominated by wiggles which modify the 
background shape. Such wiggles can easily introduce artifacts. 

In Table [T] we have summarized the S/N ratio of the three 
simulated sine waves. The S/N ratio has been defined as the ratio 
of the peak amplitude to the local background. It does not take 
into account any effects of stochastic noise in time series. The 
increase in the S/N after applying the BDf is clearly shown. For 
the data filtered by smoothing, the variations of S/N ratios are a 
direct consequence of the wiggles, and are very sensitive to the 
size of the boxcar which has been chosen. 



Table 1. S/N ratio in PSD of the three simulated sine waves at 
frequencies vi= 10"^, V2- 2 x 10"* and V3 = 10"^ Hz. 



Applied filter 




S/N 








V2 


^3 


No (raw data) 


4 


10 


140 


BDf (corrected) 


314 


5254 


3 X 10« 


Smoothing 


165 


268 


98 
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Fig. 3. (a) Power spectrum density of the raw time series; 
(b) PSD of the three simulated sine waves without the low- 
frequency trend; (c) PSD after a direct (uncorrected) application 
of the BDf; (d) the same spectrum corrected by the Q{cl>) factor; 
(e) PSD of time series filtered by removing a smoothing. The 
horizontal dotted lines correspond to the power of the simulated 
sine waves. 

4. Correction to the phase 

We have seen that the BDf provides the right frequencies of the 
oscillatory signals and the right amplitudes after correcting them 
using the Qioj) factor However, the BDf also modifies the phase 
of the oscillations with a phase-shift depending on oj. If we de- 
note (/) the real phase and i/' the phase observed in 5f, we can 
deduce a relationship between the two variables from equation 
®: 



LjAt + 71 71 I CJ 

2 2 \ (jJr 



(7) 



The phase correction varies linearly from -7r/2 to -7r, over the 
whole frequency range, from to the cut-off frequency, ojc- 

The modification of the phase induces a modification of the 
shape of the total signal. Figure |4] illustrates this phenomenon. 
The first plot shows the beating of two oscillations with close 
frequencies (0.07 and 0.08 Hz), that can be seen as a modu- 
lated carrier wave of frequency 0.075 Hz. In the filtered series, 
the amplitude is, of course, affected: it is multiplied by a factor 
2 sin(0.0757r) ^ 0.47. Furthermore the wave has a phase lead 
of about 7r/2 (if we consider a> cl>c)- The second plot shows 
what happens for a mode (vq = 0.07 Hz) and its first overtone 
(vi - 0.14 Hz). We clearly see that the resulting pattern is dif- 
ferent. This strong modification appears because the phases have 
been modified, even if the introduced phase-shift for both oscil- 
lations is almost the same {^i =s 1.8 and 1//2 ~ 2.0). 

5. Limitations in the application of the BDf 

As we have seen in the previous sections, the BDf can be a pow- 
erful detrender but it must not be overused. If there is no drift in 
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Fig. 4. Examples of time series before (solid lines) and after (dot- 
ted lines) applying the BDf: (a) The signal contains two purely 
sinusoidal oscillations with close frequencies vq = 0.07 Hz and 
vi = 0.08 Hz; (b) The signal contains a purely sinusoidal os- 
cillation with a frequency vq = 0.07 Hz and its first overtone 
(vi = 0.14 Hz). The sample rate is Af = 1 s. 



the original signal, one has to avoid applying this filter. Indeed, 
there is always a small residual background in 6f(cj) which is 
amplified when we correct the amplitudes by dividing the origi- 
nal power spectrum by the Q{a>) factor This effect is magnified 
at very low frequencies. 

When the BDf is applied to any time series - by construction 
- we lose one point of the signal (the last one) which is usually 
left at zero to maintain the same number of points as in the initial 
series. Another problem could arise when data gaps are present 
in the original time series. When the BDf is applied to data con- 
taining N gaps, we lose N additional points as the consequence 
of the edges of those gaps. Therefore, if N is small compared 
to the total number of points, then the effect can be neglected. 
However if it is not the case, it has to be taken into account and 
if the number is too big, another kind of filter should be used. 



6. Conclusion 

In this research note we have shown the advantage of removing 
low-frequency trends by applying the backwards difference fil- 
ter The background in the PSD produced by this kind of trend is 
reduced and the S/N ratio of the periodic signals are improved. 
Moreover, the transfer function of this filter, Q{cl>), has been de- 
duced in order to correct - by dividing the spectrum by this quan- 
tity - the power of the periodic signals in the PSD. The resulting 
amplitudes inferred from the PSD are the original ones. 

Filtering with a smoothing is often used, because it is a pow- 
erful detrending tool. However, it must be carefully applied and 
the choice of the boxcar size needs a preliminary study, accord- 
ing to frequencies of searched signals. The appearance of the 
wriggles must be taken into account, ideally by correcting the 
spectrum with the transfer function of this filter. 

Even though this presentation does not deal with effects of 
stochastic noise, it emphasizes the advantages of the BDf and 
explains the manner to use properly this filter in a manner which 
is very easy to implement. 
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